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Spectral properties and anomalous diffusion in the silver-mean (octonacci) quasicrystals in d = 
1,2,3 are investigated using numerical simulations of the return probability C(t) and the width 
of the wave packet w{t) for various values of the hopping strength v. In all dimensions we find 
C{t) ~ t~^ , with results suggesting a crossover from S < 1 to 6 = 1 when v is varied in d = 2, 3, 
which is compatible with the change of the spectral measure from singular continuous to absolute 
continuous; and we find w{t) ~ t^ with < /3(t;) < 1 corresponding to anomalous difi^usion. Results 
strongly suggest that I3{v) is independent of d. The scaling of the inverse participation ratio suggests 
that states remain delocalized even for very small hopping amplitude v. A study of the dynamics 
of initially localized wavepackets in large three-dimensional quasiperiodic structures furthermore 
reveals that wavepackets composed of eigenstates from an interval around the band edge diffuse 
faster than those composed of eigenstates from an interval of the band-center states: while the 
former diffuse anomalously, the latter appear to diffuse slower than any power law. 



I. INTRODUCTION 

The discovery of quasicrystalline solid-state struc- 
tures, characterized by the presence of long-range order 
but no translational symmetry-'^, which can be viewed 
as intermediates between crystalline and amorphous 
structures^ with icosahedral^ , dodecagonaP, decagonal, 
and octagonal^ orientational symmetry, as well as the 
classical wave propagation in quasiperiodic structures^, 
renewed the interest in studying quasiperiodic systems, 
in particular their transport properties^^^^. 

Quantities commonly studied to characterize the elec- 
tron dynamics in this context are the return probability 
C{t) and the mean-square displacement w{t), defined as 

C{t)^\f^dt'\^{r,^t')\\ (1) 

and 

N 

u;2(i)^^|r;-fonvI/(f„,i)|2, (2) 

n=l 

where ^{r^t) is a wave packet initially located at a site 

■^{r,t = Q)^5p^P,. (3) 

One of the characteristic properties of quasiperiodic 
systems are power-law asymptotic behaviors of these two 
quantities in the limit i — s- oo, when 

C{t)^t-\ (4) 

and 

w{t) ~ tl^. (5) 



The exponents satisfy < /3, (5 < 1, where the two equal- 
ities correspond to the limiting cases of, respectively, the 
absence of diffusion and the ballistic motion. Inbetween 
these two cases anomalous diffusion takes place (the clas- 
sical diffusion is a special case (3 — 1/2), leading to the 
anomalous transport via Einstein's relation for the zero- 
frequency conductivity: a ~ 2e^npD(T)., where D{t) = 
w'^{t)/t is the diffusion constant, and Eq. (5) leads to 
the generalized Drude formula a ~ 2e^nF^T^''~^, with 
T being a characteristic time beyond which the propaga- 
tion becomes diffusive due to the scattering*'^". Further- 
more, the shape of the diffusion front of the wavepacket 
has asymptotically the form of a stretched exponential 
exp \x/w{t)\~^ with 7 = 1/(1 — /3) in the limit t — s- 00^"'^. 

The exponent S characterizes the spectral properties 
and equals the correlation dimension of the spectral mea- 
sure (the local density of states) Furthermore, the 
wavepacket dynamics exhibits multiscaling, where differ- 
ent powers of displacement in Eq. (2) scale with different 
exponents /3 [16-22]. Pure-point, singular continuous and 
absolute continuous spectra imply, respectively, i5 = 0, 
< 5 < 1 and (5 = 1, but the reverse is not necessarily 
true in general: S > 0, P > characterize the continuous 
spectra, while 6 = usually means a pure-point spectrum 
and the absence of diffusion (/3 = 0), but with nontriv- 
ial exceptions'^. One-dimensional quasiperiodic systems 
commonly have singular continuous spectra^^'^*"'*, and 
may furthermore exhibit transitions from pure-point to 
absolute continuous behavior (being singular continuous 
at the transition point) as a parameter in the Hamil- 
tonian is varied, for instance in Harper's model of an 
electron in a magnetic field'^, or the kicked rotator'^*'. 

Eigenstates of quasiperiodic systems exhibit multifrac- 
tality, when the electron wavefunction is neither localized 
nor extended but exhibits self-similarity^^'^^'^^''^^ (defini- 
tions and further discussion are given in Sec. VI), similar 
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to the Anderson model of localization for disordered con- 
ductors near the localization-delocalization transition^^. 

The Hamiltonian H studied in this work (defined 
in Sec. II) describes a system of atoms of a d- 
dimensional product of mth approximants of the silver- 
mean (octonacci) off-diagonal model. 

In a previous work^^, the direct sum of one-dimensional 
quasiperiodic systems was studied, and it was shown 
that already in this case C{t) displays non-trivial 
dimensionally-dependent behavior. On the other hand, 
it is easy to prove analytically exactly that in quasiperi- 
odic systems defined as a direct sum of one-dimensional 
quasiperiodic systems, w{t) does not depend on the di- 
mensionality of such tilings. In the product systems stud- 
ied here this is not the case anymore and therefore they 
can be viewed as a natural next choice for the construc- 
tion of interesting higher-dimensional quasiperiodic sys- 
tems in the study of w{t). 

The rest of the paper is organized as follows: In Sec. II 
we define the systems studied, whose eigenproblem is dis- 
cussed in Sec. Ill; in Sees. IV and V results of the calcu- 
lation of C{t) and w{t) are presented and discussed, re- 
spectively; in Sec. VI participation ratios of eigenstates 
are analyzed, and finally, conclusions are presented in 
Sec. VII. 



II. DEFINITIONS 




FIG. 1: Labyrinth tiling corresponding to the Hamiltonian (8) 
for m = 5,d = 2. The distances L and S of the silver-mean 
chain are taken as 1 and v = 0.5. Then the bond lengths of 
the labyrinth tiling are \/2, v'T+l^, and v\/2. The bipartite 
structure is indicated by the open and closed circles for the 
sites, taking into account that hopping occurs only among the 
sites of different kind. The second, equivalent, labyrinth tiling 
is shown in the figure by dots. 



The silver-mean chain is defined over an alphabet 
{L, S} by the inflation rule 



P ■■ 



L 



LSL 
L 



(6) 



iterated m times starting with letter 5, with the corre- 
sponding Hamiltonian 



{nn') 



(7) 



where the sum is restricted to nearest neighbors and the 
hopping integrals t„n> take values 1 and < ii < 1 for, 
respectively, letters L and S ("large" and "small" hop- 
pings) of the letter sequence of the approximant (with 
open boundary conditions). Among various number- 
theoretical properties of the approximants, we men- 
tion that the mth approximant 's length Nm satisfies 
limm^oo Nm+i/Nm = 1 + \/2 = 5$, the so-called silver 
mean. 

The d-dimensional quasiperiodic tilings we study are 
derived from the direct product of silver-mean chains^^'^^ 
with the Hamiltonian 



(8) 



This has a particularly simple geometrical interpreta- 
tion, since it describes a particle on a rf-dimensional 
cube, with coordinates r = {xi, . . . , x^}, hopping only 



along the main diagonals of the cube: {xi, . . . ,Xd} 
{xi ± 1, . . . , ± 1}. Since the nearest neighbors of 
the cube cannot be connected by any number of such 
hops from Hm\ the Hamiltonian (8) decomposes into 
N = 2^^"^ parts defined on corresponding interpenetrat- 
ing disconnected quasiperiodic tilings. It can be shown 
that all of these Hamiltonian parts are equivalent, and 
that the tilings they define are symmetry related, due to 
the reflection symmetry of the silver-mean chain about 
the middle bond of the chain. 

Such tilings in d = 2 are known as labyrinth tilings, and 
Figure 1 shows one example, an m = 5 tiling. The tiling 
has bipartite structure. Furthermore, a two-dimensional 
projection of the silver-mean tiling in d = 3, for instance, 
looks exactly as the tiling in Figure 1, because the whole 
tiling has a layered structure: all symbols of the same 
kind in the figure belong to one layer. The distance be- 
tween successive layers are, of course, again determined 
by the silver-mean sequence. The number of sites N in 
various tilings is given in the Table I. 



III. EIGENPROBLEM 

The eigenproblem of is 

ff^lV'i) =^i|V'») , i = l,...,Ar„, (9) 
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m 


d= 1 


d = 2 


d = 3 


2 


4 


8 


16 


3 


8 


32 


128 


4 


18 


162 


1458 


5 


42 


882 


18522 


6 


100 


5000 


250000 


7 


240 


28800 


3456000 


8 


578 


167042 


48275138 


9 


1394 


971618 


677217746 


10 


3364 


5658248 


9517173136 


11 


8120 


32967200 


133846832000 


12 


19602 


192119202 


1882960298802 



TABLE I: The number of sites A'^ of the d-dimensional tihngs 
of various approximants m of the silver- mean (octonacci) quar- 
sicrystal. 

to which the eigenproblem of Hm^ reduces: 

ij = 0,...,Nm. (10) 

Regarding the subtiUng eigenproblem, its eigenstates 

are obtained simply by projecting the eigenstates of 
Hm^ onto the subtiling (with normalization constant 
N~^/'^) due to the equivalence of the individual eigen- 
problems (9). Thus, if P denotes the projector onto 
the subspace of one subtiling, the corresponding sub- 
tiling Hamiltonian is PHm^P and its eigenstates are 
_/V~i/2p|-0j^) (g) . . . (g) lipf^) with the same eigencnergics as 
in Eq. (10). The reduction of the problem to quasipcri- 
odic subtilings of the product tiling therefore only re- 
duces the degeneracy of individual eigenstates of Hm\ 
while the set of distinct eigenenergies is exactly the same. 
Furthermore, since all the Hamiltonians studied here 
have only off-diagonal matrix elements, their eigenener- 
gies necessarily come in opposite-energy pairs, and thus 
we can restrict our discussion to the spectral properties 
for states with £■ < 0. 

Figure 2 shows the dependence of the eigenenergies on 
the parameter v. A characteristic feature of the spectra 
is the existence of gaps of various sizes. In d = 1, the gap- 
labeling theorem gives an enumeration of all the possible 
gaps in the spectrum^'^''^^. For d > 1 these gaps seem to 
close for a sufficiently large v, and Fig. 2 suggests that 
gaps persist only for v < 0.6, 0.4 in d = 2, 3, respectively. 

For d > 1 the gaps close for smaller values of v due 
to the level crossings that are a consequence of the fac- 
torizability of the eigenenergies in Eq. (10). Due to this 
we expect that the spectrum will acquire finite Lebesgue 
measure and change from from fractal to continuous for 
intermediate values of v, similar to the results of Ref. [27]. 

A more difficult question, that we address in the next 
section, is whether the spectrum is singular or absolute 
continuous or a mixture of the two in different regions of 
the spectrum for various values of v. To this end we first 



note that 2 states corresponding to E = —1 for f = 
are a consequence of the open boundary conditions used 
in this work, corresponding to the electron localized, for 
small V, at the ends of the chain (in d = 1 and corners 
of the cube for d > 1), and spreading across the system 
as f ^ 1. In this respect the systems studied here are 
different from those of Ref. [35] , where two E = states 
have been induced by adding two sites at each end of the 
chain with tp ~ 0. However, for the results presented 
here this difference in boundary conditions does not play 
any important role. 

IV. RETURN PROBABILITY AND SCALING 
OF THE SPECTRAL MEASURE 

In this section we study the scaling of the return prob- 
ability C{t) by evaluating numerically Eq. (1) for the 
eigenstates of a silver-mean tiling in d = 1,2,3 con- 
structed, as described in the previous section, from eigen- 
states of the chain, obtained from the numerical diagonal- 
ization of the Hm, and taking into account that the time 
evolution is given by V-'j(i) = '/'j cxj){—iEjt). We present 
the results in Fig. 3. For small values of v there are 
pronounced steps, and this limit we analyze elsewhere^^. 

For each of the studied cases in the figure, there are 
three characteristic time-scales of the particle dynamics: 

(i) short times, characterized by a high probability to find 
the particle at the initial site, which corresponds to the 
regime when the wavepacket has only begun spreading; 

(ii) intermediate times, when C{t) does seem to behave 
according to the power law (4); and (iii) the long-time 
limit when there is a crossover into a constant C{t), due 
to the finite spatial extent of the studied systems. 

The quantitative determination of 6 is rather difficult, 
due to the large subdominant contributions to the asymp- 
totic power-law behavior^^. These appear through the 
dependence of 6 on the interval of t values taken for 
the least mean-squares linear fit in the log-log plot of 
C{t). Lacking rigorous results, one attempt to charac- 
terize these corrections to Eq. (4) would be to consider a 
class of power-law subdominant corrections by assuming 

C{t) r^t-^ + cit-^' , 0<S<S'<1. (11) 

It was noted in Ref. [24] that, in the case of the Fibonacci 

model, C{t) can be accurately fitted to Eq. (11) imder the 
assumption that S' is exactly equal to 1 , and the authors 
gave some plausible but non-rigorous arguments for this 
assumption. On the other hand, numerical verification 
of Eq. (11) by a least-squares nonlinear fit requires that 
C{t) is calculated for many decades of Int. In addition, 
verification of (5' « 1 by such a procedure would be par- 
ticularly difficult near the transition where the fit would 
be expected to yield already S ^ 1, i.e. to decompose 
the fitted function into two nearly identical terms. Fi- 
nally, the parameter v would be expected to play some 
role also in the subdominant terms, for instance through 
the dependence of the value of 6' on v. 



4 



Figiare2.jpg 



FIG. 2: Dependence of eigenenergies on the parameter d in d = 1, 2, 3, from the bottom panel to the top, respectively, for the 
iterants m = 10, 7, 6. Only eigenenergies E < are shown due to the symmetry of the spectra about the band center. 



An additional difficulty was pointed out in Ref. [36], 
where the scaling exponent of the second moment of the 
spectral measure was considered (which is known^^ to be 
equal to 6) for several quasiperiodic systems, including 
Fibonacci chains, and an apparently irregular behavior of 
the moment was found in all of the non-Fibonacci models. 
This led the authors to use a fitting procedure that does 
not involve fitting of a straight line and even to consider 
the possibility that the exponent S might not exist in 
these non-Fibonacci systems. 

In this work we fit the power law (4) to the calculated 
values of C{ti) (for many values of ti equidistant in the 
logarithmic scale) by fitting of the expression 

C(t) = 74i-'^+«'"*, (12) 

in all intervals containing ti values spanning at least one 

decade in t (for easier comparison, this corresponds to a 
change in Int of at least 2.3 in the plots), and selecting 



among the obtained fits of 6 those that have k values 
closest to 0. Should several distinct values of S occur 
that all have small k, it seems reasonable that those val- 
ues corresponding to later times are closer to the correct 
asymptotic value of S, just as it is reasonable to expect 
that for larger t values Eq. (11) becomes close to the 
power law (4). Our fitting procedure does not assume 
anything in particular about the subdominant terms of 
the powerlaw asymptotic behavior, but rather estimates 
S by requiring the absence of certain types of subdomi- 
nant terms, as quantified by k. 

We show such obtained fits in Fig. 4 for the case 
m = 6, d = 3. This is the most difficult system among 

the three cases presented in Fig. 3 because the duration 
of the intermediate time regime (ii) , where the power law 
asymptotic behavior should be most likely expected, de- 
pends only on the linear size of the system, determined 
by TO, and not on the dimensionality d. The obtained val- 
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FIG. 3: Calculated values of the return probability C{t) in d = 1, 2, and 3, from left panel to the right, respectively, for 
V = 0.1, 0.2, . . . , 0.9, respectively, from the top curve to the bottom. The dashed line in each panel has the slope —1. 




FIG. 4: Return probability in d = 3 (points) together with fit- 
ted lines close to the power law from the procedure described 
in the text. 



ues of K in the fits presented in Fig. 4 range from 10 ^ 
to 10-^. 

A characteristic feature of the fits is that for v > 0.5, 



the time intervals where the power law is obeyed best 
(i.e. where k is smallest) are shorter than the intervals 
for V < 0.4, even though for large v the curves for In C{t) 
do not show oscillatory behavior and appear to be more 
similar to the straight lines. There is, however, a small 
curvature in the dependence of \nC{t) on Int for the 
crossover times inbctwccn regimes (i) and (ii) as well as 
(ii) and (iii). This is hard to see in the plots, since it 
involves large relative changes in small curvatures, but 
easy to distinguish numerically in our fitting procedure. 

The results for S obtained from the above described 
fitting procedure applied to the results from Fig. 3 are 
presented in Fig. 5. The plot for d = 2 shows S{v) also 
for m = 6, which is systematically larger than 6{v) of the 
m = 8 system for v < Vc and smaller for v > Vc, from 
which we speculate on the possible form of the limiting 
curve of the infinite system as indicated in the figure for 
d= 2, 3. This extrapolation o{6{v) 'md= 2, 3 is in agree- 
ment with semi-quantitative considerations of this limit 
in Ref. [24] . For d= 1 the results suggest that 5 < 1 for 
all V < 1 also in the infinite system. The exponents thus 
suggest that there is a singular continuous spectrum for 
all V values in d = 1 and a possible transition from singu- 
lar continuous to absolute continuous spectra at Vc ~ 0.6 
in d = 2 (in agreement with a previous estimate'^^), and 
at Vc « 0.4-0.5 in d = 3. 

A difficulty in obtaining exact values of Vc could be 
due to the possibility that, for v !v Vc, the spectrum 
is a mixture of singular continuous and absolute continu- 
ous parts, when absolute continuous parts appear in only 
small intervals of the full spectrum that grow larger with 
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FIG. 5: Calculated values for the exponent S characterizing the asymptotic behavior C(t) ~ t"^ . Error bars are of the size of 
the symbols. Connecting lines are included for guiding the eye. The possible behavior of S in the limit of large m is indicated 
by dotted lines (see text for further discussion). 



increasing v, and when the S values obtained for v away 
from Vc are closer to the correct asymptotic values. 

A previous study^^ of the propagation of the projec- 
tion of the wave packet (3) onto various segments of 
the spectrum of the silver-mean model in d = 2 found 
that C{t) of such restricted wave packets still obeyed the 
power law, with exponents differing only slightly from 
S obtained when the full spectrum is used. We were 
able to extend this result by proving the following: If 
the spectrum is divided into a finite number of M seg- 
ments such that 1 = where Pa is the projector 
onto the subspace corresponding to the segment a, and if 
Ca{t) ~ t~^'' in each of the M segments, where Ca{t) is 
defined as in Eq. (1) with the wavepacket (3) replaced by 
its projection P„^'(r, 0), and if Eq. (4) also holds, then 
S > maxa{Sa}- 

V. ANOMALOUS DIFFUSION 

The main advantage of calculating the width w{t) of 
a wave packet is that it is more directly related to the 
transport properties than C{t), as discussed in the intro- 
duction, and that its asymptotic power-law behavior (5) 
seems to be less influenced by subdominant contribu- 
tions, and therefore easier to determine from numerical 
studies of finite-size systems, as we discuss below. 

One of the main difficulties in numerical studies of the 
anomalous diffusion is that in order to obtain w{t) ac- 
curately, many eigenvectors are needed. Calculating its 
behavior at large times also requires large system sizes. 
This puts significant constraints on the system sizes that 
can be investigated numerically, and several approxima- 
tions have been developed to circumvent these difficul- 
ties. One of them is to expand the evolution operator in 
terms of small-time increments^", or to study the time 



evolution of the position operator expanded in terms of 
Chebyshev polynomials combined with energy filtering 
via a Gaussian operator centered at a given energy and 
of the width of, e.g., 1% of the total bandwidth'*-'^ . In the 
latter work, generalized quasiperiodic Rauzi tilings have 
been studied and the authors were interested how the 
topological connectivity of the tiling influences transport 
properties of the quasicrystal in d = 2,3. The simplify- 
ing feature of the model is that the Hamiltonian is sparse 
with no matrix elements different from 1, which signifi- 
cantly speeds up any algorithm where multiplication with 
the Hamiltonian is a time-consuming part of the calcu- 
lation, like the algorithm of [41] as well as, for instance, 
the Lanczos algorithm. Another approach was to study 
the evolution of a wave packet '^a{r,t = 0) constructed 
from a subset of eigenstates from an interval a of the full 
spectrum^^, 

*„(r,t = 0) aPaJ^.,^-. (13) 

Using such wave packets, the largest system studied in 
this work is m = 7 for ti = 3 with 3456000 sites. 

Figure 6 presents results of the calculations of the 
spreading width iu{t) of a wave packet initially localized 
on a single site using all the eigenstates of the spectrum. 
For small v values pronounced steps appear, related to 
the similar steps in Fig. 3 and the structure of very nar- 
row bands of the spectra (cp. Fig. 2), that we discuss 
elsewhere^^. 

The three regimes of wave-packet dynamics discussed 
in the previous section here correspond to (i) the nearly 
ballistic propagation for early times t <1, (ii) anomalous 
diffusion, in intermediate intervals of t, characterized by 
Eq. (2) and (iii) a stationary regime of approximately 
constant w{t), due to the finite extent of the studied sys- 
tems, for large t. The maximum of w which would be 
possible for a system of linear size L can be estimated 
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FIG. 6: Calculated values of the width w{t) of the wave packet for three systems as indicated in each panel, for v = 
0.9, 0.8, . . . , 0.1, going from top curve to the bottom, respectively. The dashed line in each panel has slope 1 correspond- 
ing to the ballistic motion. 



as ^/dL/2, corresponding to a state completely local- 
ized at the corners of the sample; for the system sizes in 
Fig. 6, maximal values of Inw are approx. 7.4, 5.1, 4.5 in 
d = 1, 2, 3, respectively. On the other hand, a completely 
delocalized state would lead to Inw ~ 6.878, 4.587, 3.919, 
respectively; thus Fig. 6 shows that for large values of 
V the wave packet becomes smeared out over the entire 
sample in the regime (ill). The crossover from regime (ii) 
to (iii) becomes longer for larger d, which is due to the 
more complicated reflection of the wave packet off the 
boundaries in higher-dimensional tilings. 

The values of the exponent /3, which were obtained 
using the same fitting procedure as in the previous sec- 
tion, are presented in Fig. 7. They seem to be inde- 
pendent of the dimensionality of the system. Since the 
calculation of w{t) is computationally slower than that 
of C{t), system sizes considered here are smaller than in 
the previous section, and in particular, for v = 0.9, the 
power law was observed for only about a half a decade 
of time for d = 2, m = 7, as opposed to two decades for 
d = l,m — 10, and we could not obtain a reliable expo- 
nent using our fitting procedure for d = i,m = 6,v = 0.9. 
In the fitting procedures used in the previous section, for 
comparison, the calculated exponents were obeyed for at 
least a decade (and up to 6 decades in some cases for 
smaller v values). 

In order to assess the role of topology of the quasiperi- 
odic tiling in d = 3, we make a comparison of the wave- 
packet dynamics studied here with those of Ref. [41]. 
Since in the latter work energy-filtering was used as 
described above, we study the spreading of the wave 
packet (3) composed only of eigenstates with energy E 
in the interval < — \Ea\ < AEa, where AEa is 
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FIG. 7: Values of the exponent /9 governing anomalous diffu- 
sion of the wave packet, as extracted from the data in Fig. 6. 
The error bars are of the size of the symbols. 



chosen such that the interval contains about 2% of the 
total number of states. 

Figure 8 presents calculated values for several ener- 
gies Ea between the band edge and the band center. 
Apparently when the wave packet is composed only of 
band-edge states, the diffusion is anomalous with the ex- 
ponent very close to the value /3 ~ 0.66 obtained when 
the full spectrum is considered, as indicated in the fig- 
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lire. For other energies, there seems to be an intermediate 
time regime when anomalous diffusion also takes place, 
with the same or smaller exponents as compared to the 
one found at the band edge. Furthermore, a third kind 
of behavior also occurs, most notably at the band cen- 
ter, where the wave packet seems to spread slower than 
any power law. This is in qualitative agreement with the 
findings of [41] in the sense that there wave packets made 
out of states from the band edge also spread faster but, in 
contradistinction with the finding here, ballistic motion 
was found at the band edge and anomalous diffusion at 
the band center. 

It is unclear to us at present what kind of eigenstates 
are responsible for the slow diffusion seen at the band 
center. Qualitatively, marginally critical states^^ might 
be perhaps related to this kind of wavepacket dynamics. 

Such constructed wavepackets will inevitably have 
some "ripples" that depend on how well eigenstates from 
the chosen interval a approximate ^f.Fg, which in turn 
depends on the choice of the starting site tq. This can 
be characterized by po = \Pa6f,ra\'^, and states in Fig. 8 
have po ~ 14% and 0.7% at the band edge and the band 
center, respectively. In order to show some evidence that 
the obtained dynamics in the two cases does not depend 
on the choice of the initial site, we present w{t) in Fig. 9 
when the initial site is chosen two sites along the main 
diagonal of the system away from the center, in which 
case po « 0.3% and 11% at the band edge and the band 
center, respectively. The dependence of the values of po 
on the choice of the initial site is reflected in the range of 
values that w{t) takes, but does not seem to affect the two 
different types of dynamics. Although the contribution 
Po of the band edge states to the wave packet in Fig. 9 
is much smaller than to that in Fig. 8, the spreading of 
the entire wave packet is again characterized by the slope 
(3 Ri 0.66 in agreement with the value shown in Fig. 7. 

VI. PARTICIPATION RATIOS 

The results of the previous section strongly suggest 
that (3{v) is independent of d for the class of systems 
Hm^ studied here. In this section we attempt to link 
that with the eigenstate properties by the analysis of the 
participation number. 

The (generalized) inverse participation numbers, 
Zg{ip), are defined as 

Z,(V') = ^|V(r)r«. (14) 

r 

The participation number P2(V') = ^2^{'^)^ fo^^ instance, 
characterizes on how many sites a given state V is sig- 
nificantly different from zero: for a particle completely 
localized at a single site P2{_5r,fa) = Ij while for a Bloch 
wave P2(exp(ifc • r)/\/N) = N. The participation ratio, 
P2{ip) = P2{'4')/N, gives the fraction of the total number 
of sites where ^jJ is significantly different from 0. 



In the case of the eigenstates of Hm , we have 

ZgQiJi,) ^ ■ ■ ■ ^ = ^•••^|V'ii(a^i)---^i.(a;d)|'« 

Xl Xd 

= Zq{'4>,^)---Zq{lp,^), (15) 

or, in other words, Zq of the product state is equal to 
the product of the inverse participation numbers of the 
one-dimensional states. 

On the other hand, the inequality 

{niinZq{iji)f < Zg{i,i,) ■ ■ ■ Zgiij,,) < (maxZ,(V'i))^ 

(16) 

gives us insight into the nature of eigenstates of Hm"^: 
they are always less or equally extended than the most 
extended state of the product, and more or equally ex- 
tended than the least extended state of the product. In 
particular, they can be anisotropically extended, which 
happens whenever states from the product have quite dif- 
ferent participation numbers. For each such state of -ffm , 
there is another one with the same energy and anisotropy 
in a different direction due to the product structure of the 
Hamiltonian. 

A previous numerical study"^^ indicates that the scal- 
ing exponents of Z2 of the ground state and the band 
center state are only slightly different. Here, we investi- 
gate the scaling of the participation ratio p averaged over 
the whole spectrum as a function of Nm and v in d= 1, 
and the results are presented in Fig. 10. 

The results indicate that the average participation ra- 
tio scales with a power law in the number of sites of the 
chain, 

Pm{v) ~ Nliv), (17) 

and the inset of Fig. 10 shows the values of ^{v) obtained 
from the least mean-squares linear fit to the data for the 
largest system sizes in the figure. The results suggest that 
7(z; 1) = 0, which is not surprising since in this limit 
the quasiperiodic tiling becomes a periodic chain (we note 
that in our case p(v — !■ 1) < 1 because we have real 
instead of complex eigenstates so that the Bloch waves 
have nodal structure). On the other end, however, j(v — » 
0) = — 0.52±0.03, even though for i> = the chain breaks 
up into clusters of atoms, and therefore ^{v = 0) = —1. 
This particular limit we address elscwhere^^. For d > 1 
we have also calculated 7(1;) numerically for several v 
values and it remains the same within the error bars, 
which we attribute to the multiplicative nature of the 
participation numbers as refiected in Eqs. (15) and (16). 
This can be related to the independence of (3{v) with 
respect to the dimensionality determined in the previous 
section. 
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FIG. 8: Calculated values of the width w{t) of the initially localized wave packet composed only of states with eigenenergies 
from small energy intervals containing around 2% of the total number of states, for d = 3, m = 7, w = 0.5, (with 3456000 sites) 
going from the band edges (top curve) to the band center (bottom curve). Individual curves are shifted vertically for clarity. 
The dashed line near the topmost curve has the slope /? m 0.66 that was found in the m = 6 system using the whole spectrum. 




FIG. 9: Same as Fig. 8 but for a slightly different initial site of the wave packet (see text for further details). The dashed line 
has the slope f3 « 0.66, same as in Fig. 8. 
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!; = 0.975 




Number of sites 

FIG. 10: Dependence of the average participation ratio 
P2{v) on the number of sites Nm of the silver-mean chain 
for m = 4, . . . , 12 in a doubly-logarithmic plot for v = 
0.075, 0.1, 0.125, . . . , 0.975. The inset shows the values of the 
exponent f{v) obtained from the fit of p ~ N"'. 

VII. CONCLUSIONS 

We studied spectral properties, dynamics of wave pack- 
ets and scaling properties of cigcnstates in 1-, 2- and 3- 
dimensional systems obtained as direct products of the 
silver-mean chains, by investigating the scaling exponents 
S, P and 7 describing the asymptotic properties of, respec- 
tively, the return probability, the spreading of the wave 
packet, and the average participation ratio of eigenstates. 

The obtained values of S are compatible with the spec- 
tral measure being singular continuous in one dimension 
and undergoing a transition from singular to absolute 
continuous in d = 2,3 with large subdominant contri- 
butions. The latter would appear (if such a transition 



indeed occurs) as a systematic shift of the value of S{v) 
when the system size increases. The comparison of the 
results for S{v) for m = 6, and 8 in the two-dimensional 
systems (cp. Fig. 5) shows such a shift of the order of 
up to 10% near the speculated transition point Vc ~ 0.6 
and thus corroborates such conclusion. A more quanti- 
tative characterization of these corrections, such as, for 
instance, whether they are of the form (11), is beyond 
the scope of this work. The obtained values for f3 arc 
independent of the dimensionality, which we have linked 
to the properties of the inverse participation numbers. 

Even though these results are certainly related to the 
product structure of the quasiperiodic systems studied 
here, it is quite possible that they describe some features 
of generic higher-dimensional quasiperiodic systems. In 
particular, we found an exact relation among the expo- 
nent S and corresponding exponents Sa of the dynamics 
of the projection onto subintervals of the spectrum of a 
wave packet initially localized at a single site, and com- 
pared the dynamics of wave packets in the silver-mean 
quasiperiodic tiling constructed from about 2% of the 
total number of states near various energies in the band 
with a study of generalized quasiperiodic Rauzi tilings in 
Ref. [41]. There seems to be a similarity between the two 
in so far that wave packets of states near the band edge 
are spreading much faster than those made out of the 
states near the band center. However, while in the lat- 
ter work dynamics of these wave packets is, respectively, 
ballistic and anomalously diffusive, we find that the dy- 
namics is, respectively, anomalously diffusive and, for the 
wave packets made out of band-center states, slower than 
any anomalous diffusion. 
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